Study of maximal bipartite entanglement in valence-bond states 
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We study maximal bipartite entanglement in valence-bond (VB) states and show that the average 
bipartite entanglement E%, between a sub-system of two spins and the rest of the system, can 
be maximized through a homogenized superposition of the VB states. Our derived maximal E 2 
rapidly increases with system size and saturates at its maximum allowed value. We adopt two 
ways of generating maximal E„ states: (i) starting with a general superposition of VB states, we 
impose homogeneity; (ii) by considering a general homogeneous state, we generate isotropy. By 
using these two approaches, we construct explicitly four-qubit and six-qubit highly entangled states 
that are both isotropic and homogeneous. We also demonstrate that our maximal E% states are 
ground states of an isotropic infinite range Heisenberg model (IIRHM) and represent a new class of 
resonating-valence-bond (RVB) states. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.67.Bg, 75.10.Jm 



I. INTRODUCTION 

Quantum entanglement, a manifestation of non- 
locality, is a precious resource for quantum computation 
and quantum information [l[ and signifies correlations 
in many-body systems. Quantum algorithms that would 
significantly accelerate a classical computation must rely 
on highly entangled states since slightly entangled states 
can be simulated efficiently on a classical computer Q. 
The strength of correlations of fluctuations of observ- 
ables (such as density, magnetization, etc.) in a many- 
body system is a reflection of the degree of entanglement 
(for pure states) p|. Thus, characterization of multi- 
particle entanglement and production of maximal/high 
multi-qubit entanglement is vital for quantum compu- 
tational studies and for mutual enrichment of quantum 
information and many-body condensed matter physics. 

Intensive work on entanglement during the past decade 
has led to the proposal of numerous measures of entan- 
glement [3, @- While two-party entanglement is quite 
well understood, entanglement in a multi-party system 
is an area of immense current interest [6l4l3|. In the 
quest for maximally entangled states, so far only few- 
qubit maximally entangled states such as two-qubit Bell 
states, three-qubit Greenberger-Horne-Zeilinger (GHZ) 
states, and four-qubit Higuchi-Sudbcry (HS) states [l4[ 
have been clearly identified. It would be of considerable 
interest to generate these maximally entangled states as 
the ground states of a physically realizable spin model. 
While the GHZ states could be obtained as the ground 
state of an anisotropic Heisenberg model [l5j maximally 
entangled four-qubit and five-qubit states could be ob- 
tained only as excited states [la ]. 

VB states are useful in understanding the entangle- 
ment properties of cluster states that enable universal 
quantum computation [13, Gil- Furthermore, VB states 
are also useful as variational states with valence bond 
length being the variational degree of freedom [l9l - l2lj |. 



RVB states are the linear superpositions of VB states 
and have provided interesting insights for understand- 
ing strongly correlated phenomena such as spin liquid 
physics in frustrated magnets [22], physics of high T c 
cuprates [23U25j , superconductivity in organic solids [2(| , 
insulator-superconductor transition in boron-doped dia- 
mond [27]], etc. Furthermore, RVB states have also been 
proposed as robust basis states for topological quantum 
computation [28j . The multi- particle entanglement has 
also been investigated in RVB states that were proposed 
as states close in energy to the ground state of the Heisen- 
berg Hamiltonian [29M3l| . However, there has been no 
explicit construction of RVB states that would represent 
maximally entangled VB states. 

The main results of this paper are threefold. First, 
we study the entanglement defined by von Neumann en- 
tropy and show that isotropy and homogeneity maximize 
the bipartite entanglement E%. Second, based on these 
principles of isotropy and homogeneity, we propose two 
approaches to construct entangled states that maximize 
E%. Lastly, we analyze these states in terms of ground 
states of IIRHM and observe that these states represent 
a new class of RVB states. 

The remainder of this paper is organized as follows. 
In section II, we introduce entanglement entropy and 
demonstrate explicitly that isotropy and homogeneity 
maximize entanglement E% between two spins and the 
rest of the system in the two limits of even-N-qubits, i.e., 
for N = 4 and N — > 00. For intermediate even-N-spin 
states, our proposed entangled states are only shown to 
maximize E% among VB states. In section III, we pro- 
vide two approaches to construct maximal E% entangled 
states: first approach deals with homogenizing isotropic 
states while the second one generates isotropy in homo- 
geneous states. We also demonstrate our procedures by 
constructing maximal E% states for four-qubit and six- 
qubit systems. In section IV, we show that our E% max 
states are ground states of our IIRHM while in section V, 
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we make a comparison of our E% max states with the RVB 
states and observe that our maximal E% states represent 
a new class of RVB states. Lastly, in the final section VI, 
we conclude after commenting on the physical rcalizabil- 
ity and decoherence free nature of our proposed high E 2 V 
entangled states. 



II. ENTANGLEMENT FOR TWO-QUBIT 
REDUCED DENSITY MATRIX IN ISOTROPIC 
SYSTEM. 

For a bipartite system AB in a pure state, von Neu- 
mann entropy E v measures the entanglement between 
the subsystems A and B. From the reduced density ma- 
trices pa = trsp AB and ps = ti'Ap AB obtained from the 
pure state p AB , we obtain 

E v (pa) = -tr(p A log 2 p A ) 

= -tr(p B log 2 p B ) =E v (p B ). 



Using the basis | |), | f), and S l = 
density matrix can be written as [3 
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Throughout this paper, we consider only states 1^) that 
are eigenstates of the z-component of the total spin oper- 
ator {S^ otal ) with eigenvalue Sj.; furthermore, we focus 
on only isotropic states. It then follows that, (5f) = 
and (iSj ) = leading to the single-qubit density matrix 
to be maximally mixed and thus maximizing entangle- 
ment E v (pi). On realizing that (Sf) = 0, we obtain the 
following expression for the two-qubit reduced density 
matrix 1331: 
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Here, isotropy implies 0.5(St Sj) = 0.5(S+St) = 
(SfSJ) = (SVS?) = {SfSp. Thus, the von Neumann 
entropy E v (p^ ) can be expressed as 



E v (p i3 ) = 2 



For our states, S^ otal \^ 
(5fEjS?> = 0, that is, 

/](SjS*) = 



3(l + 4(S?S?»log 2 (l+4(S?S;» 
(l-12(£fS|))log 2 (l-12(Sf^))](4) 

= which implies that 



(5) 



We will now maximize the total entanglement entropy 
Yli jjti Ev(Pij) subject to the above constraint in Eq. (|5|). 
To this end, we will employ the method of Lagrange mul- 
tipliers and define the Lagrange function A as follows: 



A=Y > E v ( Pij )-J2\ [Y^istsr 



(6) 



Then, setting 



OA 



A, = 3 log 2 



= yields 

(1 - l2{SfSQ) 
(l + 4<Sf5*» 



(7) 



which implies that the optimal {SfS*) is independent of 
j for j ^ i. Consequently, it follows from Eq. ([5]) that 
Y.i^i E v{Pij) is maximized when (SfSj) = - 
i.e., when the isotropic state is homogeneous. The av- 
erage entanglement entropy between the subsystem of 
two spins and the rest of the system (of N — 2 spins) 
E% = [1/N(N — 1)] J2i j^i Ev(Pij) has a maximum value 
given by 
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It is interesting to note that for N — > 00, the above 
expression yields E% — > 2. In fact, i? 2 approaches the 
maximum possible value 2 quite rapidly as can be seen 
from Fig. [TJ Next, we observe that for N = 4, our ex- 
pression for E"l max in Eq. ^ yields the same entangle- 
ment entropy value 1 + 0.5 log 2 3 as that obtained for the 
four-qubit maximally entangled HS state [TFj, HH ■ Fur- 
thermore, our approach explains why all the pairs in the 
HS state give the same entanglement value. Contrast- 
ingly, for the isotropic ground state, in the case of a 
Heisenberg chain with nearest-neighbor interaction, for 
four spins the entanglement entropy E 2 = 1.21 (with 
(S*Sj) = —0.5/3) while for an infinite chain the von Neu- 
mann entropy E 2 V = 1.37 (with (SfSj) « -0.443/3) both 
of which are far less than our maximal i? 2 values above. 
It is of interest to note that, when N = 4 or N — > 00, 
the maximal values of E% for isotropic systems are the 
same as the maximal i? 2 values for general systems (i.e., 
systems that can be either isotropic or non-isotropic). 

We also note that homogeneity, under the constraint 
of Eq. (0), maximizes the i-concurrence (I c ) [3~i| given 
by 



N(N — 1) 



£</2[l-Tr(^)] 



(9) 



As shown in Fig. 1, I c also monotonically increases with 
system size. Although we considered von Neumann en- 
tropy and i-concurrence as entanglement measures, our 
homogeneous states should also maximize other measures 
of entanglement for in VB systems. 
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FIG. 1: Normalized entanglement E 2 /E^ ax , for two-qubit re- 
duced density matrix, measured by (a) von Neumann entropy 
(Ev/Ev, max ) and (b) i-concurrence (I c /Ic,max), for N-qubit VB 
systems. 



III. MAXIMAL Et ENTANGLED STATES 

In this section we offer two different approaches for 
constructing maximal E% states; the first approach in- 
volves producing homogeneity in isotropic states while 
the second deals with generating isotropy in homoge- 
neous states. As examples of our prescribed procedure, 
we construct highly entangled states for four and six 
qubits. 



A. Generating homogeneity in isotropic states 

In this section we will proceed to construct entangled 
states for N-spins that maximize E%. We first note that 
there are (N — 1)!! states with St = which are the 
product of N/2 two-spin singlet states of the form | 'hj- 
) — I it) (with no pair of singlets sharing a spin). Of these 
(JV— 1)!! product combinations with St = 0, only n Cn_ — 

n Cn_ 1 = N\/[{N/2)\{N/2 + 1)!] products are linearly 
independent. A particular set of linearly independent 
St = states are the Rumer states [HI [3(| that are 
made up of non-crossing singlets. 

Next, we demonstrate how highly entangled states can 
be constructed by starting from isotropic states and mak- 
ing them homogeneous. Since we are dealing with spin- 
singlets, we consider t spins and y \. spins in our basis 
states \tpk) (= | of 02 — 0jv) where of =t or |). Using su- 
perposition of these basis states, with all states having 
equal probability, we construct the states 
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(10) 



which wc will now prove to be homogeneous. Throughout 
this paper, for convenience, we ignore the normalization 
constants in our spin states. Let the spin at site i be ei- 
ther j" or J,; then the probability that the spin at site j ^ i 



FIG. 2: Linearly independent non-crossing Rumer diagrams 
for a four qubit system . 



is in the same state is [^/2)-\ ] +n /2 wnue the probability 



that it is in the opposite state is [^/2) 1 -i ] +n/2 ' 
fore, wc can write 



Thcrc- 



(SfSj 



— - l 

2 x 



'K - l) + £L 

2 L > ' 2 
1 



N 
2 
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2 L ) ~ 2 
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4(N-1) 



(11) 



Thus, when all the basis states are equally probable, we 
get (S?Sj) = — 4(^-1) ■ Now, we are in a position to 
form (from isotropic St = states) homogenized states 
that are highly entangled. 

Four Qubit maximal E% states: For N = 4 
spins, we have two linearly independent non-crossing 
Rumer diagrams which are shown in Fig. [2] The 
states corresponding to these diagrams can be expressed 
as |$f 2 12= °) ® l^ff = °) and |f i4 4= °) ® l$23 3 = °) whcr c 
|ff/=°) = [\ t)«| I) j - I t>,;] is a two-spin singlet 
state for spins at sites i and j with Sij being the total 
spin of Si and Sj. It is worth noting that 
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I* 
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Sl4=0\ 

14 / 



S 23 =0\ 
23 / 



This means that the crossing singlet |f fij 3 °) 
depends linearly on the two non-crossing 
] ) ® |ff 4 34=0 ) and 



If 



Sl4 = 0\ 



i*fr =0 > 



(12) 

singlets 
Using 



the above relation one can establish that there are only 
n Cn — n Cn_ 1 linearly independent St = states. 

Next, we take the following linear superposition that 
will lead us to the desired entangled state: 



|f 4 ) = r x e l 



r 2 e' 



Sl2=0\ 

12 / 

Sl4=0\ 

14 / 



i*fr =0 ) 



if 



(13) 



where r^e 1 ^ represents a general coefficient. 

Evaluating Eq. ([13]) (and setting |f 4 ) = \^ N )hom for 
N = 4), we get the various coefficients (shown in Table 
IJ) for each of the basis state \ip%) occurring in Eq. (jXOJ) . 

The expression for |f 4 ) [in Eq. (|T3|) ] assumes the 
homogenized superposed form shown in Eq. (fT0| for 
two linearly independent sets of solutions of the coeffi- 
cients in Table |U the two sets of solutions are {rie 101 = 
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Basis states \ipl) 


Coefficients e iS " in Eq. (flu). 


\nn) 


rie iSl - r 2 e i92 


\mv 


ne ifll - r 2 e i02 


imt) 


—rie 


\mi) 


-ne Wl 


\nu) 


r 2 e i02 


lutr) 


r 2 e 102 



TABLE I: Basis states \tpl) of Eq. (fTU|) for a four-qubit system 
and the corresponding coefficients e* fc obtained by setting 
|* 4 ) [in Eq. CS])] equal to |*jv)w [in Eq. {TO)] with TV = 4. 



e i2u-/3 r2 gi6i2 _ e i47r/3| anc i ^ s complex conjugate. The 
states corresponding to these solutions are 

|* 4 ) = e t2,r / 3 (|$f^ =0 )® |$ff= )) 
+ e i47r / 3 (|$f^ =0 >® |$f| 3=0 )) 

= [(i nn> + 1 mt» + ^(i mt> + 1 im» 
+«! (|ttu> + iutt))] 

= |*hs) (14) 

and its complex conjugate. In the above equation, W3 = 
e i2ir/3^ j e ^ a cu h e roo ^ f un ity. The above state is the 
well known Higuchi-Sudbcry (HS) state for four qubits 

Six Qubit maximal states: The linearly inde- 
pendent states in this case are five in number (i.e., the 
five non-crossing Rumer diagrams of Fig. |3]). Thus as in 
the four-qubit case, we begin with the following general 



Basis states tpf. 


Coefficients e iSk in Eq. ifTu). 


\nnn) 


r ie i01 -r 2 e if>2 -r 3 e i0a -r±e i0i -r 5 e i05 
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-ne' 8 ' +r A e 10 * 


ItUttl) 


-ne Wl +r 5 e i05 


imm> 


ne iSl -r 4 e iBi 


mmt) 


-ne i01 +r 4 e iBi 


mm) 


ne iBl -r 5 e if>5 


imtit) 


ne iBl -r 3 e i<>3 


\mm) 


-ne iBl +r 2 e iB2 +r 3 e lS3 +r i e lBi +r 5 e i05 


ittuti) 


+r 2 e 102 +r 4 e w * 


limit) 




\mm) 




\mm) 


-r 2 e 102 -r±e 10i 


\mm) 


—rr,e 5 


httw> 


+r 5 e i05 


imm> 


-r 2 e 102 -r 5 e i05 


limn) 


+r 2 e t02 +r 5 e i<>5 


\inm) 


-r a e i03 


imm) 


+r 3 e i03 


iimtt) 


-r 2 e 102 -rs,e i03 


\nnu-) 


+r 2 e 102 +r a e i03 



TABLE II: Basis states of Eq. (fTOjl for a six-qubit system 
and the corresponding coefficients e lSk determined from the 
equation |* 6 > = |* 6 ) hom . 
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(b) 

FIG. 3: Homogenized linear combination of the five non- 
crossing St — singlet states (or non-crossing Rumer dia- 
grams) in (a) and (b) give maximal E% entanglement for six 
qubits. The state \^c) m Eq. (|19[) is a linear combination of 
the three diagrams in (a); on the other hand, |$®) and |$®) 
in Eqs. (|16p and (|17|1 . respectively, represent different linear 
combinations of all the three diagrams in (b). 



linear superposition: 



= ne* 




=°>« 


5 |<£ff 4=0 ) c 


3 \t>!r°) 


+ r 2 e l ° 2 
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5 i*fr =0 ) 



and find solutions that homogenize \^f 6 ), i.e., make |\l/ 6 ) 
assume the \^e}hom form in Eq. (fT0|) . Next, we evaluate 
Eq. (fT5|) and get the various coefficients for each basis 
state mentioned in Eq. (fTU|) . The expressions for 
the coefficients e lSk in Eq. (JTUJ) arc shown in Table [TTJ 
As in the four-qubit case, we make the state in Eq. (jT5]) 
homogeneous, i.e., assume the form given in Eq. (fTUl) . 

We observe that, since there are only five linearly inde- 
pendent Rumer states, we can construct at most five lin- 
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early independent superpositions of the form |\1/ 6 ) given 
by Eq. (|15[) ; thus we can expect five independent sets 
of coefficients r^e 4 ^ producing homogeneity. Here, we 
present five linearly independent solutions: 

1. {n = r 2 = 0;(r 3 e' 93 ,r 4 e^,r 5 e w ") = (w^wf)}, 

2. {n = r 2 = 0;(r 3 e' 9 \r 4 e* e \r 5 e^) = (wj.wf.wf)}, 

3. {n =r 2 =0;(r 3 e ld3 ,r 4 e l9 \r 5 e ld ") = (-1,+1,-1)}, 

4. { ri e m = -1 + e M ; r 2 e ?;e2 = f + e M ; 

(r 3 e^,r 4 e^,r 5 e^) = (-f,-f,-f)}, 



5. {ne 101 = -1 



f e-'Vae'* 2 = 1 + e" 1 "; 
(r 3 e^,r 4 e^,r 5 e^) = (-l,-f,-f)}. 



where w 4 = e l27r / 4 is a fourth root of unity and e IQ is 
a complex number. In arriving at the above solutions, 
use has been made of the fact that when Ylj=i &u3ti = 0, 
we have two pairs of unit vectors with each pair giving 
a zero sum (as shown in Fig. 2] for one the three pos- 
sibilities). The states corresponding to the first and the 
second solutions are 
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S 56 =0 
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)) 



+w2(|*S' =0 ) ® l*2? = °) ® l$fj 4= °», (16) 
and its complex conjugate, while the state that corre- 



sponds to the third solution is 

|*g) = -(|<&? 
+(\* 
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The states that are obtained from the fourth and the 
fifth solutions can be written as 
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® l*25 B 


=°) 


® i* it 


= 0) - 





and its complex conjugate. Using e la 
state in Eq. (fT8|) can be rewritten as 



-<(= 



(18) 
-0J4,), the 



I* 



^ 4 (|$f 2 12=0 ) € 
+^l(l^fi 4=0 ) 

■+^2(i*lr=°) 



8 l*2 S | 5=0 ) 



1$ 



I* 



23 



65 

i$fr=°)). (19) 



=0 » 



This state represents the superposition of the three 
Rumcr diagrams of Fig. [UJa) with the second term in 
the sum representing crossing singlets. 

Evaluating the tensor products in the above Eqs. (|16p . 
(fl~7|) . and (fl"9|) . we rewrite these states as 



\*i) = [(I mm) - 1 mm)) +m\ mm) + 1 mm) + 1 imti) - 1 mm) - 1 mm) - \ mm)) 
+ loK\ mm) + 1 mm) + 1 utm) - 1 mm) - 1 mm) - 1 umt» 
+ mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm))], 



(20) 



I*?) 



and 



\**) = 



mm) - 1 mm)) - (i mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) 
+ (i mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 wmt» 
- (i mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm))], 



mm)) 



(21) 



mm) - 1 mm)) + mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 
+oj 2 4 (\ mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 
+ojK\ mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm))}- 



(22) 



B. Producing isotropy in homogeneous states 



value Ef, 



The condition of isotropy in a homoge- 



Next, we produce isotropy in a given homogeneous 
state; such states should also yield the entanglement 
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neous state can be expressed as Yli S" |*iv)/K>ro = with 
a = x 7 y,z; this implies that |*7v)/iom = 0. Now, 

5 —0 S- —0 

reflection about z-axis for produces —^ij 3 ', then 

in an isotropic St = state comprising of singlets, re- 
flection leads to a coefficient of (— 1) a" for the parent 
state. Hence, in the homogeneous state of Eq. (fT0|). im- 
posing isotropy implies that flipping all spins should pro- 
duce a coefficient of (— l)" 5 " for the parent state; in fact, 
the coefficients e lSk of the basis state \a\ <7f...erj^) and 
its spin flipped version differ only by a factor (— 
Thus the number of unknown coefficients e lSk in Eq. 

(fTU|) is reduced to 0.5 ^Csj; this also implies that 

12i S~\^N)hom = yields the same equations as the con- 
dition S^\^N)hom = 0. Therefore, we can write the 
four-qubit homogeneous state as 



+ e*[imt) + um)] 
+ e*[itm) + mtt)]. 



(23) 



Furthermore, the condition of isotropy in the above 
four-qubit homogeneous state yields 



( e #i + e i<t>2 + e «fa) [| + | 



; = 1 



= 0. 



This leads to the expression e 1 



mt> + 1 ittt>] 

(24) 



J4>3 



0, i.e., 

zero valued coefficients for all the four basis states. The 
only two possible solutions to this equation (up to a phase 
factor) are {e** 1 = l,e^ 2 = w 3 ,e^ 3 = w§} and its com- 
plex conjugate. The states obtained from these solutions 
are the same as those obtained earlier by imposing homo- 
geneity on isotropic states and represent the four-qubit 
HS-state and its complex conjugate. 

Next, we consider the case of six-qubit homogeneous 
state which can be written as: 



I* 



6 / horn 



e i*i [ 


mm) - 


mm)] 


e *<fe [ 


mm) - 


mm)} 


e i<t>3 [ 


mm) - 


mm)] 


g^4 [ 


mm) - 


mm)} 


e#5 [ 


mm) - 


umt>] 




mm) - 


utm>] 




mm) - 


mm)] 


e i4>s [ 


mm) - 


mm)} 




mm) - 


nun)] 




[| mitt) - 


I mm)} 



(25) 



The isotropy condition for the above six-qubit homoge- 
neous state (expressed as $t \^ n) hom = 0) yields 
the various basis states and the corresponding coefficients 
shown in table IIIII The set of fifteen equations (cor- 



Basis states 

dJllLdllltJU 111 


Zero valued coefficients of 

Uclolo a LcLLco U U LcllllfciLl 11 (Jill 

using isotropy on * 6 ) tora . 


itttm) 


g«0I + + e '4>s _ e i4>r 


lt;tm> 


e '<Pl _|_ g*<fe _ gJ06 _ g*<£lO 


Itmtt) 


g*01 _ g* 1 ^! _ g 4< ^8 _ p i <i>9 


imm> 


g'^2 _|_ g*>>4 J_ g»09 _ g*<^>6 


limn) 


g»<#>2 _ g*<#>3 _ gI08 1 g*01O 


imm> 


e i<A 3 + gi<#>4 + gi05 + g^8 


\mm) 


_ e iH + e i<?4 + e i4>7 + e i<Pw 


IttUtt) 


g^5 + g^6 _ e i4>9 _|_ g^lO 


limit) 


_g'^3 _ g»04 1 g!06 _ g*07 


imtit) 


_g^l _ gi02 + gi06 + g^S 


imttt) 


_gi</>2 + e i4>3 + e i4>7 _ R i<t>% 


\urm) 


_g*</>5 _ gi'fre _|_ g'^7 _ g*<#>8 


itttm> 


-e i<A7 + e i0s + e^ 9 - e^ 10 


litmt) 


_gi<^l _ g»03 _ gi05 _J_ g*<#>9 


Ittmt) 


— g'^2 _ g»04 _ gi<fe _ g*<#>10 



TABLE III: The basis states and the corresponding zero 
valued coefficients obtained from the condition of isotropy 

Etl S 4 + |*e>fcom=0. 



responding to the zero valued coefficients of the basis 
states) in table IIIII can be reduced to the following set of 
five linearly independent equations: 



e ^ 



e i<f>z 

c*04 







„ J <?9 



g'^io = o 



0. 



(26) 
(27) 
(28) 
(29) 
(30) 



The above equations can be viewed as the sum of unit 
vectors on a circle in the complex plane. Since the sum 



is of the form J2j=i e% 



0, it implies the conditions 



shown in Fig. 01 The five linearly independent solutions 
of the above set of Eqs. ([26])-(|30]) are given in Table HVl 
This set of solutions produces the same entangled states 
iM), I*6>. I*c>. and (!*«))*] as obtained in the 

previous section. 



C. Discussion and comparison 

We will now discuss generating E% max RVB states for 
the general case of an even N spin system. The num- 
ber of coefficients rie %Si , needed to generate a super- 
posed state |* w ) using all non-crossing Rumer states, 
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FIG. 4: Geometric depiction of zero sum of four unit vectors 
on a unit circle in the complex plane. This represents one of 
the three possibilities of the zero sum £j =1 e* a J = 0. The 



other two possibilities are e %ai + e"* 2 
and e iai + e iai = 0, e iQ2 + e lQ3 = 0. 



0, e' Q3 + e" 



0: 



is the same as the total number of linearly independent 
Rumcr states n Cn_ — iv C»_ 1 . The number of indcpcn- 

2 2 1 

dent basis states \ip%) in Eq. (HOJ is 0.5 (^Cf)- O n 

setting 1^^) = \^N)hom and equating the coefficients of 

the various basis states, we get 0.5 f^Cwj equations; 

from these equations, on eliminating r^e* 61 ' in terms of 

the various e lSk . we get —0.5 {^Cn \ + n Cn__ 1 number 

of linearly independent equations [of the form of Eqs. 

([26 |) -(|30 ]) when N = 6] in terms of 0.5 f w C|j number 

of unknown coefficients e k . Thus, we expect the num- 
ber of independent E% max RVB states to be equal to 



N C n — N C n (i.e., the number of unknown e l5k minus 

2 2 1 x 

the number of independent equations) which is the same 
as the number of linearly independent Rumer states. For 
instance, for N = 8, we get 21 linearly independent equa- 
tions in terms of 35 different e lSk ; thus the number of 
independent E% max RVB states is 14 (= 35 — 21) which 
is the total number of non-crossing Rumer states. 

Next, we observe that the von Neumann entropy 
for two-qubit reduced density matrix obtained from 
(l*a))*> l*b), l*c>) or (l*c))*i is 1-921964 which is the 
same as E% max proposed by our general formula in Eq. 
((8]). The conjecture made by Brown et al. [H[ that the 
multi-qubit maximally entangled states always have their 
reduced single qubit density matrix maximally mixed is 
satisfied by our states \^ 6 ). It should be noted that the 
highly entangled six-qubit state reported in Ref. [l(| , al- 
though yields a higher entanglement, is not an eigenstatc 
of the S^nt r i operator. 

Cabello [331 nas a ^ so constructed supersinglets of four 
and six qubits which are decoherence free. However, 
the average two particle von Neumann entanglement en- 
tropy E% for these supersinglets is less than the entan- 
glement for our states constructed above. In the case 

(2)\ 

of the four qubit supersinglet state 04 > , the entan- 
glement entropies for all possible bipartitions are given 
by E v ( Pl2 ) = E V ( P34 ) = 1.5849, E v ( Pl3 ) = E v (p 2A ) = 
1.2075, E v (p 14 ) = E v (p 23 ) = 1.2075 which together yield 
the average two particle entropy value 1.3333, i.e., a 
quantity clearly smaller than the von Neumann entropy 
El = 1 + 0.5 log 2 3 « 1.79248 obtained for the four-qubit 
maximally entangled HS state. In the six qubit supersin- 
(2) \ 

glet SI' 3 ) case, the various bipartitions produce the av- 
erage entropy E% = 1.657996 which is noticeably smaller 
than the entropy value E% max = 1.921964 obtained for 
our six qubit entangled states. 



Coefficients in 
Eqs. |26]l-P0)l 


Five linearly independent solutions 


Sol. 1 


Sol. 2 


Sol. 3 


Sol. 4 


Sol. 5 


e i<h 


1 


1 


1 


1 


1 


e i4>i 






-1 


-e ia 


-e~ ia 


e i<t>3 


w| 




-1 


-e la 


-e~ ia 


e i04 


2 

— LU4 


*2 
— CJ 4 


-1 


e ia 




e i05 


2 




1 


e la 


e -ia 




2 

— (xJ 4 


*2 
— CJ4 


-1 


1 


1 


e^ 7 




-wf 


1 


1 


1 








1 


-e ia 


-e~ ia 




— U>4 


-U)l 


1 


1 


1 




— LUi 


— U]% 


1 


-e ia 


-e~ ia 



TABLE IV: Sets of five linearly independent solutions of the 
Eqs. (|2oT)-PJl). 



IV. GENERATING HIGHLY ENTANGLED 
GROUND STATES USING IIRHM 

In this section we will demonstrate that the E% max 
entangled states obtained in the previous section from 
the homogenization of isotropic states (or from imposing 
isotropy on homogeneous states) are the ground states 
of a spin Hamiltonian. Thus we begin with the IIRHM 
Hamiltonian 



-Hiirhm 



i,j>i 



, (31) 



and show that certain superpositions of the ground states 
of IIRHM will produce the same amount of entanglement 
as given by Eq. ((8j) . We observe that [S^ otal , -ffiiRHM 1 



and [S% otal , Hubum] = 0. In Eq. 



Totals 

H, we take J = 
J* /{N — 1) (where J* is a finite quantity) so that the 
energy per site remains finite as N — >• 00. The eigenstates 
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of i?iiRHM correspond to eigenenergies 



E 



St 



S T {S T + 1) 



37V 



(32) 



where St is the total spin eigenvalue. The ground state 
corresponds to St = which is rotationally invariant and 
also implies that S T = 0. Next, we note the interesting 
fact that 



[Si. S3 + S2.S4 + S1.S4 + S2.S3 



1$ 



S 3 4.=0\ 

34 / 



0. 



using which it can be shown that any St — state made 
up of N/ 2 singlets will be an eigenstate of Eq. (|3T|) . Thus 
the Rumer states are ground states of J?iir,hm • Hence, for 
a homogenized linear combination of non-crossing Rumer 
states, the entanglement is given by Eq. (JSJ. 

It is important to note that, while all possible St = 
states are ground states of IIRHM, for some Hamiltoni- 
ans only some of the possible St = states are ground 
states. For instance, for the nearest-neighbor antiferro- 
magnetic Heisenberg Hamiltonian or for the Majumdar- 
Ghosh model Hamiltonian [38| only some of the possible 
St = states are ground states and consequently (for 
these systems) it is not possible to construct a homog- 
enized linear combination of St = ground states that 
produces maximal possible entanglement E% max given 
byEq. m. 



V. RESONATING- VALENCE-BOND PICTURE 

The construction of our entangled maximal E* states 
that form the ground states of IIRHM can be visualized 
by using a RVB picture. Our maximal E% states can be 
regarded as a new class of RVB states made of homog- 
enized superposition of isotropic St = valence bond 
states. We will now compare the entanglement proper- 
ties of our RVB states and the general RVB states j^rvb 
of Ref. given below: 



vb 



E 

i a eA 



f(ii, ...,i M ,ji, —,jM)\(ii,ji)—(iM,3M)), 



where M represents the number of sites in each sub- 
lattice and / is assumed to be isotropic over the lattice. 
Also, |(i fc ,jfc)} = ^(| t>ij - I i)i k \ t)iJ denotes 
the singlet dimer connecting a site in sub-lattice A with 
a site in sub-lattice B. The valence bond basis (used for 
the above RVB state |\P)rvb) form an over complete set 
while our RVB states are constructed from a complete 



set of N C 



N 



Ci 



states. 



The rotational invariance of the two qubit reduced den- 
sity matrix of the RVB states allows us to write them in 
the form of a Werner state: 



Pw{p) =p\(ij)){(ij)\ 



1 



P 



h 



(33) 



where for 1/3 < p < 1 the Werner state has the spins 
at i and j entangled with each other. For the special 
case of the "RVB gas" for the l^rvb state, where / is 
a constant (corresponding to equal amplitude superpo- 
sition of all bipartite valence bond coverings), one gets 
the exact result p = 1/3 + 2/{3N) which implies that 
all finite size systems have a non-zero tangle (or en- 
tanglement) between the two sites [3(3]. Next, for the 
"RVB liquid" case (involving equal amplitude superpo- 
sition of all nearest-neighbor singlet valence bond cover- 
ings of a lattice) , Monte Carlo calculations extrapolated 
to the thermodynamic limit yield p = 0.3946(3) > 1/3 
[3f3 |. i.e., a non-zero tangle between the two sites [33 |. 
In contrast, our maximal E% RVB states yield zero en- 
tanglement between the two spins for all system sizes as 
demonstrated below. It has been shown that the SU(2) 
symmetry of the RVB states ensures that the two-spin 
correlation function and the parameter p of the Werner 
state are related as 



(*|5 < .S J -|*) = --p. 
Then, since our RVB states produce 

/ 01 qx\ / qV qV\ / qz qz 



1 



4(N - I)' 



(34) 



(35) 



it follows that p = jtzt an d thus for all even N > 4 
systems we get zero ent ang lement between the two sites. 
Lastly, based on Ref. [30, we find that the monogamy 
argument yields the bound p < 1/3 + 2/ — 1) 
while the quantum telecloning argument produces the 
bound p < 1/3 + 2/[3(N — 1)]. Compared to our exact 
value of p = 1/(N — 1), both these bounds are weaker 
bounds (i.e., show zero two-site entanglement with cer- 
tainty only in the thermodynamic limit). Thus (among 
various RVB states) we see that our entangled high E% 
RVB states, while producing maximum entanglement be- 
tween a pair and the rest of the system, yield zero entan- 
glement among the two spins of the pair. 

We will now remark on the entanglement of any set of 
n (> 2) sites with the rest of the system. The valence 
bond entanglement entropy picture, developed in refer- 
ences (HJ and [HI , uses the number of valence bond two- 
spin singlets shared by the two subsystems as a measure 
of their entanglement with each other. Thus, our 'ho- 
mogenized' RVB states, will always have shared bonds 
between the two subsystems and thus show high bipar- 
tite entanglement. However, exact quantitative analysis 
needs to be carried out. Extending our approach to de- 
riving the n-qubit reduced density matrix in terms of n- 
particle correlation functions (when n > 2) and obtaining 
entanglement expressions and wavefunctions correspond- 
ing to maximal bipartite entanglement is a non-trivial 
exercise and is left for future studies. In this connection, 
our notion of homogeneity [as expressed in Eq. (ITU1) ] is 
useful even for higher order correlation functions (since 
\^N)hom generates correlation functions (Sf S| ...Sf S^} 
that are independent of the site indices i,j,...,l,m when 
no two site indices respresent the same site). 
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VI. CONCLUSIONS. 

Although both spins in a two-spin singlet state | 
) — I 4-t) are monogamous (i.e., they cannot be entangled 
with any other spin), by using a homogenized superposi- 
tion of valence bond states (each of which is a product of 
N/2 two-spin singlets), we managed to distribute entan- 
glement efficiently such that any pair is maximally en- 
tangled with the rest of the RVB system while concomi- 
tantly the constituent spins of the pair are completely 
unentangled with each other. Thus we get the converse 
of monogamy for a pair of spins! Now, while total spin 
zero states are quite commonly ground states (as shown 
by Lieb-Mattis theorem it has not been recognized 

that one can generate high bipartite E% entanglement 
from such states. Our maximal E% RVB states are phys- 
ically realizable in systems such as infinite-range large 
U jt Hubbard model and infinite-range hard-core boson 



model with frustrated hopping [4lJ when they are at half- 
filling. 

It can be shown that the eigenstates of our Hamilto- 
nian IIRHM are decoherence free when the system is cou- 
pled to a bath of either local optical phonons or global 
optical phonons [42|. Since our RVB states (that are 
constructed from homogenization of VB states) form the 
ground states of IIRHM, they are decoherence free. 

In summary, our maximal E% RVB states have led us to 
IIRHM for their realization and IIRHM in turn generated 
maximal E% states as highly robust ground states. 
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